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Vortex ring instability and its sound 


By R. Verzicco 1 AND K. Shariff 2 


This work carries earlier finite-difference calculations of the Widnall instability of 
vortex rings into the late non-linear stage. Plots of energy in azimuthal Fourier 
modes indicate that low-order modes dominate at large times; their structure and 
dynamics remain unexplored, however. An attempt was made to calculate the 
acoustic signal using the theory of Mohring (1978), valid for unbounded flow. This 
theory shows that only low-order azimuthal modes contribute to the sound. As 
a check on the effects of axial periodicity and a slip wall at large radius imposed 
by the numerical scheme, the acoustic integrals were also computed in a truncated 
region. Half of the terms contributing to the sound have large differences between 
the two regions, and the results are therefore unreliable. The error is less severe for 
a contribution involving only the m = 2 mode, and its low frequency is consistent 
with a free elliptic bending wave on a thin ring. 


1. Introduction 

Many shear flows of practical interest contain large-scale coherent motions, and 
some aspects of their dynamics must be responsible for the generation of sound 
(smaller scale motions affect the dynamics of the large scales but themselves radi- 
ate as higher order poles). What these aspects might be is difficult to identify from 
experiments since the most energetic or most visually apparent motions may not be 
the most radiative and a detailed knowledge of time- dependence is required. It is 
a usual procedure in fluid dynamics to extract from a flow the most representative 
structure with the hope that its evolution alone still retains the important features 
of the whole flow. In the present case the evolution of azimuthal instabilities on an 
isolated vortex ring is numerically simulated into the late non-linear stage. We are 
motivated by the following: (i) Observations using phased arrays of microphones 
which indicate that the acoustic source location in jets is near the end of the po- 
tential core where vortex ring structures break down (Bridges & Hussain 1987, p. 
309). Implicit in the consideration of a single ring is the premise that the actual 
breakdown process does not involve interactions between rings. Factors which may 
be important in reality are enhanced growth-rate and mode selection due to mu- 
tual straining and tearing or pairing, (ii) A good fit (Shariff & Leonard 1992, p. 
275) of the acoustic frequency of a turbulent vortex ring measured by Zaitsev et 
al. (1990) to the axisymmetric elliptic core ring but a poor fit to the elliptic bending 
mode. Further studies along these lines were carried out in Zaitsev & Ko’pev (1993) 
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FIGURE 1. Acoustic signal due to axisymmetric elliptically distorted vortex ring. 
, T/u = 1200; , T/u = 3000. 

and Kop’ev & Chernyshev (1993). (iii) The observations of Maxworthy (1977) of 
the development of swirling flow and a solitary bulge wave following the Widnall 
instability. 


2* Methods 

Calculations were performed using a second order finite difference code in cylin- 
drical coordinates (x,r, </>) assuming axial periodicity (with length L x ) and a slip 
wall at r — R 0 as described in Verzicco & Orlandi (1993). The spurious straining 
due to the infinite row of rings and due to the wall is estimated to be less than 2% 
of the strain (due to ring curvature) which drives the instability. Units are chosen 
so that the circulation, P, and ring radius, i?, are unity. To avoid effects of axisym- 
metric unsteadiness, a ring with an initially Gaussian core (of radius a = .4131) 
and T/u — 3000 was axisymmetrically relaxed to a quasi steady state in which 
becomes some function of the streamfunction. The number of grid intervals in each 
direction is 128 with a smaller radial spacing, obtained by a non uniform mesh, 
near the vortex core. The domain size is R 0 = L x = 6. A random divergence- free 
perturbation is applied to the basic-state velocity. 

The acoustic pressure, p 0 , is obtained from the theory of Mohring (1978): 

Pa = |§ - A/c o), Qij = ^ J *.(* X u)j d 3 x , (1) 

where A = XiXi . The triple prime denotes the third time derivative and is numer- 
ically evaluated using cubic splines. Due to the factor X{Xj only the symmetrized 
Qij are required, and in cylindrical coordinates one obtains: 
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t = 0 t = 15 t = 25 

FIGURE 2. Contour- plots of azimuthal vorticity fon an elliptically distorted Gaus- 
sian core vortex ring at T/u = 3000. Contour increments are A = 0.1, , axis 

of symmetry. The ring motion is downward. 


Qn = 2nli (0 mode only) 

Q22 = ~ (A'i - h ~ 2/i — is) (0 and 2 modes) 

Q33 = ^ (A - 2/i + I& - A'i) (0 and 2 modes) 

Q12 + Q21 = * (h -J1-J2 + K 3 ) (1 mode only) 
Q13 + Q31 = it ( Ji - J 3 - h + K 2 ) (1 mode only) 
Q23 + Q32 = - h - A' 5 ) (2 mode only) 


The various integrals in Eq. (2) are 


71 -j 

1 xr 2 £ 4 ( 0 ) dxdr, 

ii 

f xr 2 Q x (l)dxdr, I 3 = j 

[ xr 2 uj x (l)dxdr , 

II 

1 xr 2 Q r (2)dxdr, 

ii 

f xr 2 u> r { 2) dxdr, I 2 = J 

^ xr 2 u5^>(2) dxdr , 

II 

1 xr 2 u; ( i > (2)dxdr , 

Ji = 

* 

j x 2 ri2 r (l)dxdr, J 2 = 

1 x 2 ruj<f>(l)dxdr , 

J 3 = 

J 

1 x 2 rQ^(l)dxdr, 

^ = 

j x 2 ru r (l)dxdr, K\ = 

J r 3 u) x (2) dxdr, 

II 

. cs 

j r z u^{\)dxdr, 

ii 

£ 

f r 3 u> ,*,(1) dxdr, K 3 = J 

p 

r 3 u5 x (2) dxdr , 


where, for instance, Q z {m) and w I (m) denote the mth sin and cos azimuthal modes, 
respectively, of the axial vorticity. Note that only the modes m < 2 are involved. 
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Figure 3. Evolution of energy in azimuthal modes. Mode index m: • , 0; ■ , 1; 
o , 2; □ , 3; ^ , 4; o , 5; v , 6; + , 7. (a) Total energy. ( b ) Energy in azimuthal 
velocity component. Certain modes have been omitted for clarity. 

The theory of Mohring is valid for an unbounded domain in which the vorticity 
decays exponentially at infinity. This holds provided the computational domain 
is large compared with the vortical region. As the ring leaves and re-enters the 
period, what its position would be in an unbounded domain is tracked in order to 
obtain an x value in the unbounded domain for every grid point. The field is phase 
shifted axially so that the vortical region is centered in the integration domain and 
mid-point rule quadrature is employed. To assess errors due to out-lying vorticity, 
integration is also performed in a smaller domain of radial and axial dimensions 
R t = 5ii 0 /6 and Lt = 51^/6, respectively. As a test, the acoustic signal for an 
axisymmetric ring with an elliptically distorted Gaussian core was obtained (Fig. 1). 
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t = 156 * = 214 

FIGURE 4. Contours of vorticity magnitude in the ( r y 0 ) plane of maximum | u |. 
The spatial resolution of the simulation is twice the resolution of the figure. 


It has the form of a damped cosine. The damping is due to the readjustment of the 
initial vorticity distribution towards an equilibrium distribution (Fig. 2). For an 
inviscid axisymmetric elliptic core ring with uniform vorticity, the acoustic signal is 
easily worked out (Shariff et al. 1989, p. 112) from the dynamics of Moore (1980). 
For the present case it is an undamped cosine with amplitude .022 which is close to 
the initial amplitude computed. Also, the period of the acoustic signal agrees with 
the theoretical prediction. 
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FIGURE 5. Acoustic components. 
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Full-domain; , Truncated domain. 


3* Results 

Applying the inviscid stability theory of Widnall et ai (1974), valid for thin 
cores, to the Gaussian profile, one finds that the so-called ‘non- rotating second 
radial mode’ corresponds to m = 2.26R/& = 5.47 so that either m = 5 or 6 should 
be most unstable. The theory also predicts that there should be a band of growing 
modes around the most unstable one. Fig. 3a shows that the numerical result is 
consistent with this picture. Since the initial perturbation is not a combination of 
eigenmodes, there is an initial transient, after which the m = 6 mode is dominant 
showing a large range of linear growth. At t = 150 the m = 6 mode gives way to 
m — 4, likely due to viscous spreading of a (<); this switching of the dominant mode 
has also been previously observed by Hasselbrink (1992) and Shariff et al. (1994). 
At later times the non-linearly amplified m = 1 and 2 modes become dominant 
and all modes saturate. That the overall ring shape follows the dominant mode is 
shown in Fig. 4 where contour plots of |u>| in the (r, <j>) plane of maximum |u>| are 
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displayed. 

Maxworthy (1977) used dye visualization to infer that following breaking of the 
instability wave, a swirling flow developed in the azimuthal direction accompanied 
by a propagating solitary bulge wave. Fig. 36 plots modal energies in the azimuthal 
velocity component to illustrate the rapid development of a mean swirl (• ). This 
fact together with the nonlinear growth of the m = 1 wave already shown in Fig. 3a 
may be related to Maxworthy’s observation. 

Fig. 5 shows contributions to the acoustic signal during the linear instability 
phase comparing the results for the full and truncated integration domains. For 
the diagonal terms the two domains have differences comparable to the amplitude. 
These terms involve the ra = 0 and m = 2 modes. The off-diagonal terms have 
smaller differences mainly in the form of high frequency oscillations present in the 
truncated domain. The particle turn-around time is about 12 initially, and it cor- 
responds to the period of acoustic oscillations due to an elliptically distorted core. 
On the other hand, the term Q 23 + Q ^ r 2 , which involves only the m = 2 mode, has a 
larger period. This is consistent with the frequency ratio of the two modes for thin 
uniform vorticity cores: using Eq. (35) in Widnall & Sullivan (1973) we obtain 

log(a/J2), 

/elliptic core +1 

which gives an estimate of 0.3 for the present case ( a is the radius of the uniform 
vorticity core). 

Further work using the same method should either focus on the sound for axi- 
ally periodic flow (cylindrical wavefronts at infinity) using the appropriate acoustic 
analogy or employ a larger domain. 

4. Conclusions 

The time evolution of azimuthal instabilities on a thick isolated vortex ring into 
the late non-linear stage was studied by a numerical simulation of the Navier-Stokes 
equations. Initially the m = 6 mode is most amplified, consistent with the linear 
theory of Widnall et ai (1974). Subsequently, due to the viscous spreading of the 
ring, the m = 6 gave way to m = 4. In the late non-linear stage low order modes 
(m = 2 and m — 1) became dominant and a rapid growth of a mean swirl is present, 
which may be related to some dye observations of Maxworthy (1977). 

The acoustic signal computed using the theory of Mohring (1978) was presented 
for the linear instability phase. Since this theory is valid for unbounded flow, each 
contribution (Qij) was computed in a full and in a truncated domain to estimate 
errors associated with outlying vorticity. The errors were quite large for the diagonal 
terms. Further work is needed to eliminate these errors. 
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